AN OPTIMAL TRANSPORT VIEW ON SCHRODINGER'S EQUATION 

Max-K. von Renesse || 

Abstract 

We show that the Schrodinger equation is a lift of Newton's law of motion V^ /i = — V F(fj,) on 
the space of probability measures, where derivatives are taken w.r.t. the Wasserstein Riemannian 
metric. Here the potential \x — > F([i) is the sum of the total classical potential energy (V, jjl) of the 
extended system and its Fisher information 4- J* |Vln/z| 2 ci/i. The precise relation is established 
via a well known ('Madelung') transform which is shown to be a symplectic submersion of the 
<^i \ standard symplectic structure of complex valued functions into the canonical symplectic space 

over the Wasserstein space. All computations are conducted in the framework of Otto's formal 
Riemannian calculus for optimal transportation of probability measures. 

a: 

1 INTRODUCTION 

— Recent applications of optimal transport theory have demonstrated that certain analytical and geo- 
i — eaetric problems on finite dimensional Riemannian manifolds (M, g) or more general metric measure 

fiaces {X,d,m) can nicely be treated in the corresponding ('Wasserstein') space of probability measures 
2 (X) = {/j£ V{X) I j x d 2 (x, o)fi(dx) < 00} equipped with the quadratic Wasserstein metric 
^— > . 

<& • C r r ^j 1/2 

Bdw(fi, is) = mi I / / d 2 (x, y)U(dx, dy) | II G V(X 2 ), U(X x A) = u, U(A x X) = (i(A), A G B(X) \ . 

^T!his metric corresponds to a relaxed version of Monge's optimal transportation problem with cost function 
i-efay) = d 2 (x,y) 

^O '■ inf J J c(x, Ty)fi(dx) \ T : X -► X, T^ = v \ , 

32Jyith T*n denoting the image (push forward) measure of \x G V(X) under the map T. 

OQTlhe physical relevance of the Wasserstein distance was highlighted by the works of e.g. BENAMOU- 
ToRENIER [1] and OTTO [12] who established in the smooth Riemannian case X = M and smooth initial 
. distribution \x 

■ ' <(,,,) = inf \ f 1 [ \VM*)\Wdx)dt ^^(]0,l[xM),^ Ml6 C7([0,l],P(M)) 

I Jo Jm 



£ 



fit = -div(V0i/it),t G]0, l[,/i = /i,/ii = v 
showing that dy^ is associated to a formal Riemannian structure on V(M) given by 



T^V(M) = {tp : M -f R, f ip(x)dx = 0}, 

Jm 

f |V0| 2 t//i, for i) = -div(/iV0). 
Jm 
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In view of the continuity equation 

fit = - div($i/x t ) 

for a smooth flow (t,x) — > $t(x) on M, acting on measures 11 through push forward Lit = ($ t )*/io, this 
identifies the Riemannian energy of a curve t — » lit G V(M) with the minimal required kinetic energy 

E oM= / HAsllr 7?(M) ds = / / \^(x,s)\ 2 Li s (dx)ds. 
Jo s Jo Jm 

A major reason for the success of this framework is the interpretation of evolution equations of type 

d t u = div(«tVF'(V)), 

with F' being the L 2 -Frechet derivative of some smooth functional F on L 2 (M, dx), as c?w-gradient ('steep- 
est descent') flow 

fi = -V w FGu) 

for the measures ii(dx) = u(x)dx. Properties of the flow my thus be deduced from the geometry of the fun- 
tional F with respect to dyy. A particularly important case is the Boltzmann entropy F(u) = J M ulnudx 
which induces the heat flow. 

In this note we propose an example of another natural class of dynamical systems associated with the 
Riemannian metric on V(M) and which can be written as 

v^ = -v w nu). (i) 

Equation ([T]) describes the Hamiltonian flow on TV(M) induced from the Lagrangian 

L F :TP(M)^R; L F (ff>) = \ |H|^ - F(/x) for ^ 6 T„V{M) 

with the functional F : V(M) — > R now playing the role of a potential field for the system. Apart from the 
closely related recent work [TT] it seems that a systematic approach to such Hamiltonian flows on V(M) 
is missing in the literature. The example we want to propose is obtained by choosing 

Fiji) = f V(x)n{dx) + ^/(/i), (2) 

Jm o 

where 

-^(/-0 = / |Vln//| 2 (i/i. 
Jm 

We show that via an appropriate transform the flow jT]) solves the Schrodinger equation 

iHdt® = -h 2 /2A^ + V V. (3) 

The functional I is known today as Fisher information. Physically I (/a) is the instantaneous kinetic 
energy required by the unperturbed heat flow at state ii. The prominent role of / for quantum behaviour 



was noticed long ago, e.g. in a classical paper by BOHM [3J, using the following well-known system of 
generalized Hamilton- Jacobi and transport equations 

dtS + 1 -\VS\ 2 + V + y (iVln/il 2 - %) = 

dtp + div(/iV,S) = 0. 

This system was poposed by MADELUNG very early [5] as an equivalent description of the wave function 
\& = ^//lefi 5 under the Schrodinger equation. In the sequel it will be referred to as Madelung flow. Various 
attempts to derive it from first order principles can be found in the physics literature, e.g. most recently 
in0. 

Our present note starts with the observation that equations (J3J) and ([T]) are essentially the same (theorem 
12.11) . where the latter is understood in the sense of LOTT's recently proposed second order calculus on 
Wasserstein space, c.f. [8]. A virtue of formula (CQ) is its very intuitive physical interpretation as Newton's 
law for the motion of an extended system with inertia (we have put mass density equal to one) . Accelera- 
tion comes from a a gradient field of a potential F which is the total mechanical potential of the extended 
system plus its 'kinetic potential' w.r.t. the heat flow. (As usual the case of a classical single particle 
moving in a potential field is embedded naturally in ([T]) if one puts h = and /x = S x .) 

Secondly we show that the two equations (CQ) and ([3]) are, modulo constant phase shifts, symplectically 
equivalent. More precisely, we compute the canonical symplectic form on the tangent bundle TV(M) 
induced from the Levi-Civita connection of the Wasserstein metric on V(M) and show that the map 
\l/ = |\l/|eh 5 i— > — div(|\l/| 2 VS'), which we shall call Madelung transform, is a symplectic submersion of 
the standard Hamiltonian structure of the Schrodinger equation on the space of complex valued functions 
into the Hamiltonian structure associated to (CQ) on the tangent bundle TV(M). Except for its curiosity 
in Wasserstein geometry this result seems to support the point of view of some authors that the familiar 
complex valued form (J3J of the Schrodinger equation is the consequence of a smart choice of coordinates 
in which the intuitive but unhandy dynamical system ([1]) resp. (J3J) can be solved very efficiently 

Obviously, much of what is presented below resembles the familiar Schrodinger folklore, c.f. in particular 
NELSON's theory of stochastic mechanics [10] and its follow-ups, e.g. [15]. And in fact nothing really new 
about the Schrodinger equation itself is implied at this point. Our aim is the connection to Wasserstein ge- 
ometry which in our view gives a very intuitive picture. Finally, we emphasize that all of our computations 
are completely formal, a rigorous mathematical treatment of these ideas is subject to future work. 



2 SCHRODINGER EQUATION FROM NEWTON's LAW OF MOTION ON (P(M), d 



\V) 



The computations below are conducted on the formal Riemannian manifold of fully supported smooth 
probability measures equipped with the Wasserstein metric tensor, as initiated in [121 03] and extended in 
[H], ignoring full mathematical generality or rigor. (The basic background material taken from [HI H2] can 
be found in the appendix.) In the sequel we shall often identify /i G V°°(M) with its density fi = dfi/dx. 

Theorem 2.1. For V G C°°(M) let F : V°°(M) — > M. defined as in (TJ)]. Then any smooth local solution 



t — > /i(t) G V(M) of ([T]) yields a local solution (n t , S t ) of the Madelung flow (jlj), where 

S(x,t) = S(x,t) + / L F (S a , fi a )da 
Jo 

and S(x,t) is the velocity potential of the flow fi, i.e. satisfying J M Sdfi = and fi t = — div( VS^/i). 
Conversely, let (/i f , S t ) be smooth a local solution of (Jlj) then t — > /x t G V(M) solves (jTj). 

Proof. Let yU solve (jTj) where V w is the Wasserstein gradient and V^/i is the covariant derivative 

associated to the Levi-Civita connection on TV(M). Let (x,t) — > S(x,t) denote the velocity potential of 
\i (cf. section E]), then according to [HI proposition 4.24] the left hand side of (pP) is computed as 

-divLv(d t S+^\VS\ 2 

where the right hand side of (JTj) equals (cf. section [5j) 

divLv(v + y(|Vln/i| 2 --A/i)J J . 

Since \x t is fully supported on M this implies 

d t S + -\VS\ 2 + V + — (|Vln/i| 2 - -A/i) = c(t) 
2 8 /i 

for some function c(t). To compute c(t) note that due to the normalization (S^fh) — 

= dt(S t ,IH) 
^c(t)-^(\VS\ 2 ,dn)-F(^) + (S,fi) 

= <t) - \(\VS\ 2 , dp) - F(/i) + (| VS\ 2 , /i) = c(t) + L F (S t , fi t ). 

Hence the pair t — ► (St, IM) with S(x,t) = S(x,t) + J L F (S a , ix a )da solves (jlj). The converse statement is 
now also obvious. □ 

Corollary 2.2. For V G C°°(M) /ei F : V°°(M) -► K de/med as m (TJj. Taen any smooia /oca/ solution 

t -> /x(t) G 7>(M) o/ 

v v V/t = _ v w F(Ai)) 

yze/ds a /oca/ solution of the Schrodinger equation ([3]) w'a 



^(i^j^V^Mei^ (5) 

w/iere 

5(x, £) = S^x, £) + / L F (S (T ,^ a )da 
Jo 

and S(x, t) is the velocity potential of the flow fi, i.e. satisfying J M Sd/i = and fi t — — div(VS , ty u). 



Remark 2.3. The passage from S to S = S + const, does not bear any physical relevance, since two wave 
functions \1/, \1/ with \I/ = e tK ^ for some k G R parameterize the same physical system. Accordingly the 
Schrodinger equation should probably rather be understood in the sense of ihd t [^} = —h 2 /2A[^/] + [\&] 1/ 
for a flow of equivalence classes of wave functions. On the level of representatives this amounts to the 
equation 

3 «(.) : R + -»• R : z^* = -h 2 /2A$! + ^V + ik^. 

Remark 2.4. The <i>v-gradient fl° w on ^P{M) for F as in (J2J) corresponding to the overdamped limit of 
(JT]) gives a nonlinear 4th-order equation which is sometimes called the 'Derrida-Lebowitz-Speer-Spohn' or 
'quantum-drift-diffusion' equation. A rigorous treatment of it can be found in [BJ. 

The usual argument for the derivation of Euler-Lagrange equations yields the following statement. 

Corollary 2.5. ForV G C°°(M) let F : V°°(M) — > K defined as in (TJ|). Then any smooth local Lagrangian 
flow [0, e] 3 t —> fi t G TV°°(M) associated to Lp yields a local solution of the Schrodinger equation 

ihd t V = -h 2 /2A^ + <&V 

via 



V{t,x) = y/J4J~x)eW x ' t) 

where 

S(x,t) = S(x,t) + J L F (S a , fi a )da 
Jo 

and S(x, t) is the velocity potential of the flow \x, i. e. satisfying J M Sdfi = and fit — — div(VS'jyu) . 

Remark 2.6. An equivalent version of theorem l2.1l puts \& = y/]I(x, t)eTi s ^ x,t ' where t —> (— div(VS , i /i t ), fx t ) 
is a Lagrangian flow for Lp and S is chosen to satisfy for all t > 

(S t , m) - (S , Ho) = / L F (fi s )ds. 
Jo 

3 HAMILTONIAN STRUCTURE OF THE MADELUNG FLOW ON TV(M) 

In this section we show that the Madelung flow (j4]) has a Hamiltonian structure w.r.t. the canonical 
symplectic form induced from the Wasserstein metric tensor on the tangent bundle TV(M). To this aim 
we use the representation 

TV(M) = {- div(Vj» | / g C°°(M),h e V(M)}. 

Definition 3.1 (Standard Vector Fields on TV(M)). Each pair (>,</>) G C°°(M) x C°°(M) induces a 
vector field V^^ on TV{M) via 

V^(-div(V//x))=7 

where t — > ^^(t) = j(t) G TV(M) is the curve satisfying 

7 (£) = -div(/i(£)V(/ + £0)) 
Ht = exp(tV^)*/i 



Recall that the standard symplectic form on the tangent bundle of a Riemannian manifold is given by 
uj = dQ, where the canonical 1-form is defined as 

e(X) = (£,7r*(X)) T7r£ , XeT^TM), 

and where tt denotes the projection map it : TM —>■ M. 

Proposition 3.2. Let ujyv G A 2 (TV(M)) be the standard symplectic form associated to the Wasserstein 
Riemannian structure on V(M), then 

uw(V^,V^)(-div(Vf»)) = (VV, V0) M - (V^, V0) At (6) 

Proof. We use the formula 

^ W) vfo) = V^9(Vfo) - V^Q(V^) - 0([^, Vy ), (7) 

where [V^, Vj,x] denotes the Lie-bracket of the vector fields V^ and V$i. From the definition of 6 we 
obtain 

©(^)(-div(V//i)) = (V/,V^. 

Hence 

W©%*)) = 4«, e(v ^ (7 ^ (t)) 

= 4, (V(/ + t0),V^) Mt) 

at|t_o (g) 

= (V0,V^) M - / V/ • V^(- divV^/i)rfx 
= (V0,V^) At + / V(V/-V^)V^/x 

J M 

Next, since measures tangential variations only one gets that 

0([VW^D(-div(V^)) = (VJ, [V^, V$) M . (9) 

Finally, it is easy to check that 

/ V(V/-W0W>^- / V(V/-V^)V^/i-(V/,[V^,V^]) At = 0, 

which together with ((7|), (|SJ) and (|H]) establishes the claim. □ 

Remark 3.3. Proposition ^. 21 shows that cuyv is the lift of the standard symplectic form on TM to TV(M). 
This corresponds to the result in [HI section 6], which however is less explicit than formula (El). 

Using the the Riemannian inner product in each fiber of TV(M) the Hamiltonian associated with Lp is 

H F :TV(M)^R; H F (- div(V//x)) = £ / | V/|> + F(/i) (10) 

6 



Proposition 3.4. Let Xp denote the Hamiltonian vector field Xp induced on TV(M) from Hp and ujy^, 
then 

X F (-div{Vfn)) = V f _ ( i |v/|2+y+ ^ (|vlnju|2 _ 2 ^ )) (-div(V//x)) 

Proof. Fix if), cj) G C°°(M) and let V^(.) denote the corresponding standard vector field. Let t — > j(t) = 
— div((V/ + t(f>)fi t ), where fi t = exp(£V^)*//, denote the corresponding curve on TV(M), then 

V^(H F )(-div(Vfn)) = d tlt=0 H F (j(t)) 

= d tlt=0 (l [ |v(/ + ^)| 2 d^ + (F,/x t ) + ^-/(^)) 

= 1 + 11 + III, 



where 



VfVcf>dfi + l- I |V/| 2 (-div(W7i)) 



and 



m 2 j M 

(V/,V0),+ (V^V(^|V/| 2 )) 



11= f V(-div(Vipfj)) = (W,W), 

Jm 



111 = % [ 2Vln^V( div(VV,/x) )^+^ / |Vln/i| 2 (-div(V^)) 



fi 2 2A/y 

= ^«V^,V( ^)), + (V^,V|Vln/i| 2 ),) 

o /i 

Hence, collecting terms 

1 f) 2 An 

V^{H F ){- div(V//*)) = (V/, V0), - (V(-(-| V/| 2 + V + -(| V InH 2 - 2-^))), Vif)),. 

From this and formula ^ the claim follows. □ 

Corollary 3.5. T/ie pair t — >• (S t , Ht) £ C°°(M) x V{M) solves the Madelung flow equation (j4]) i/ and 
onfr/ ift—>— div(V S t fit) £ TV{M) is an integral curve for Xp. 

4 THE MADELUNG TRANSFORM AS A SYMPLECTIC SUBMERSION 

In this section we prove that the two equations (JTJ and ([3]) are related via a symplectic submersion. 

Definition 4.1. A smooth map s : (M, uj) —* (N, rf) between two symplectic manifolds is called a symplectic 
submersion if its differential s* : TM — » TN is surjective and satisfies ri(s*X,s*Y) = ou(X,Y) for all 
X,Y eTM. 



Note that this definition implies in particular that the map s itself is surjective. The following proposition 
is easily verified. Its meaning is that in order to solve a Hamiltonian system on N we may look for solutions 
for the lifted Hamiltonian g o s on the larger state space M and project them via s back again to N. 

Proposition 4.2. Let s : (M, u) — > (N, 77) be a symplectic submersion and let f £ C°°(M) and g £ C°°(N) 
with g o s = f , then s maps Hamiltonian flows associated to f on (M,ou) to Hamiltonian flows associated 
to g on (JV, 77). 

Let now C(M) = C°°(M; C) denote the linear space of smooth complex valued functions on M. Identifying 
as usual the tangent space over an element \& £ C with C, TC is naturally equipped with the symplectic 
form 



uj c (F, G) = -2 f Im(F • G)(x)dx. 

Jm 



It is a well-known fact that the Schrodinger equation (J3J) is the Hamiltonian flow induced from the sym- 
plectic form h ■ ujc and the Hamiltonian function on C 



H S (V) = — f \\7^\ 2 dx+ f \V(x)\ 2 V{x)dx. 

2 Jm Jm 



Let C*(M) denote the subset of nowhere vanishing functions from C such that j M \^(x)\ 2 dx = 1 and note 
that C*(M) is invariant under the Schrodinger flow. 

Assuming simple connectedness of M implies (via a standard lifting theorem of algebraic topology) that 
each function \l/ £ C* admits a decomposition \l/ = l^le^ 5 , where the smooth field S : M — > R is uniquely 
defined up to an additive constant h2nk, k £ N. Hence we may define a the Madelung transform 

o : C*(M) -* TP(M), a(V) = - div(|^| 2 V5). (11) 

For the next theorem recall that in our definition of TV(M) we assume that the supporting measures are 
smooth and strictly positive on M. 

Theorem 4.3. Let M be simply connected. Then the Madelung transform 

a : C*{M) -► TV(M), a(|*|e* 5 ) = - div(|*| 2 VS) 

defines symplectic submersion from (C*(M) , h ■ ujc) to (TV(M),ujy\;) which preserves the Hamiltonian, i.e. 

Hs = Hp o a. 

Remark 4.4. Together with proposition 14.21 this result presents the Schrodinger equation (jHJ) as a sym- 
plectic lifting of Newton's law on Wasserstein space (OQ) to the larger space C*(M), and which can be solved 
much easier because it is linear. Projecting the solution down to TV(M) via a yields the desired solution 
to ([1]). Going in inverse direction from ([TJ) to (J3J) requires a scalar correction term in the phase field, c.f. 
remark 12.31 

Proof of theorem^ Clearly, o"(C*(M)) = TV(M). To see that a : C*{M) -> TV(M) is a submersion fix 
a reference point £ M, then for each r £ [0, 2irh[ the map r = t^' 

t : TV(M) -> a(M), - div(VS>) -► v^e^ (s " (5(0) - r)) , 



is a bijection from TV(M) to the subset {\l/ G C*, pfr(O) = e*' - } which satisfies a o r = IcIt-p(m)- This 
proves that the differential s* of s is surjective. 

To prove that a is symplectic let ^ = ^[Jie^ G C* with /(0) = r E [0,27rh[ and let 77 = — div(/iV/) = 
ct(\I/) G TV{M). Again due to the identity a o r = Id TV ^q it suffices to prove that t*uj c = l/h ■ ujw on 
T V (TV(M)). Since the set {V^(- div(>V/)) |^,0 G C°°(M)} spans the full tangent space T V (TP(M)), 
it remains to verify 

for all ip,(j),ip,(j) G C°°(M). By definition of V^ and r = r^ for yUj := exp(tV^)*(yu) and c(t) : = 
f(0)+t<j>(0)-r 

t.Vm = dt\t=o^e i{f+t ^ cm = & (VW~ div (W/i)) + V^ ~ ' 
Hence 

wcfo^T,^) = -2 j^(-L(-div(VV>ju)) • (-^(0 + c)) 

+ ^(0 + c) • 2~p(~ div(V^//)))ete 
= ^((VV,V0) M -(V0,V^)) = ^(%V^) 
Finally, for * = r(-( divV///)), V* = ^//le^dVln/i + |V/) such that 

and J |\l/(x)| 2 l^(a;)(ia; = (V, /i) which establishes the third claim Hg = Up o a of the theorem. □ 

5 APPENDIX - BASIC FORMAL RIEMANNIAN CALCULUS ON V(M) 

Let V2{M) denote the set of Borel probability measures \x on a smooth closed finite dimensional Riemannian 
manifold (M,g) having finite second moment f M d 2 (o,x)/j,(dx) < 00. As argued in [8] the subsequent 
calculations make strict mathematical sense on the dyy-dense subset of smooth fully supported probabilities 
V°°(M) C V2{M) which shall often be identified with their corresponding density \x = dfx/dx. 

Vector Fields on V(M) and Velocity Potentials. 

A function <fi G C^°(M) induces a flow on V(M) via push forward 

t-+nt= ($ t V<?i )*Ak), 

where t — > <3>f is the local flow of difformorphisms on M induced from the vector field V</> G T(M) starting 
from $0 = Mm- The continuity equation yields the infinitesimal variation of jjl G V(M) as 

A = d t \t=oftt = - div(V0/i) G T^V). 



Hence the function induces a vector field V^ G T(V(M)) by 

^(/i) = -div(V0/i), 
acting on smooth functionals F : V(M) — > R via 

^(F)(/i) = % =0 F(/i - ediv(V#)) = % =0 F(($ 4 V V) 
with Riemannian norm 

\\V*(n)\\l v = [ \V<P\ 2 (x)fi(dx). 
Jm 

Conversely, each smooth variation if) G T^iV) can be identified with 

ip = -V^(fi) with 4> = Gpi/), 

where G^ is the Green operator for A M : — > — div(/iV0) on Lq(M, dx) = Ll(M,dx) fl {(f,dx) = 0}. 
Hence, for each ^ G T^P there exists a unique G C°° fl L 2 (M, dx) such that 

■0 = — div(/iV0) and (0, /i) = 0, 

which we call velocity potential for ip G T fJL V{M). 

Riemannian Gradient on V(M). 

The Riemannian gradient of a smooth functional F : Dom(F) C V(M) — > R is computed to be 

where x — ► DF\^{x) is the L 2 (M, dx)-Frechet-derivative of F in /i, which is defined through the relation 

d ele=0 F(fi + ef) = / DF^x)t(x)dx, 

for all £ chosen from a suitable dense set of test functions in L 2 (M, dx). The following examples are easily 

obtained. 

F(fx) = J M <}>{x)n{dx), V W F, M = VM = - div(V0/i) 

F(fx) = J M /i log fxdx, V W F| M = — div(/iVlog/i) = — A/^ 

F(/i) = f M |Vln//| 2 d/x, V W F, M = -div(/iV(|Vln/i| 2 - JA//)). 

Here A denotes the Laplace-Beltrami operator on (M,g). As a consequence, the Boltzmann entropy 
induces the heat equation as gradient flow on V(M), and the information functional is the norm-square 
of its gradient, i.e. 

||V w Ent|J^= ||-div(//Vlog/x)||^ = J |Vlog/i| 2 d/i = /(/i). 



Covariant Derivative. 

The Koszul identity for the Levi-Civita connection and a straightforward computation of commutators 

hv 



show [8] for the covariant derivative V w associtated to dw that 



(V^ 2 ,^ 3 ) Tm = / Hess 0^0!, V0 2 )d//. 

Jm 

For a smooth curve t — > ji{t) with /i 4 = V^ t this yields 
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